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Abstract A selected review of topics at the border of hard and soft 
physics is given. Particular emphasis is placed on diffraction dissociation 
at Fermilab and HERA. Recently, significant differences between diffraction 
dissociation at HERA and at Fermilab have become apparent. This may 
suggest that one already is reaching nonlinear (unitarity) effects which are 
extending from the soft physics region into the semihard regime of QCD. 



1 Introduction 

The focus in this paper is on the regime of hardness near the borderline 
between hard and soft high energy collisions with a special emphasis on 
searching for nonlinear QCD effects. This is an opportune time for such a 
discussion as there is now a significant body of complementary data from 
deep inelastic scattering and from hadron-hadron collisions. Perhaps the 
major object of this review is to compare and contrast hadronic and deep 
inelastic collisions, especially diffraction dissociation where there are major 
differences between the hadronic and virtual photon initiated processes. 

Sec. 2 is devoted to a brief review of some soft physics results on total cross 
sections and diffraction dissociation. Despite the fact that total cross sections 
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grow[l] as {s/soY, with e ~ 0.1, through the highest Fermilab energies it is 
argued that there already is strong evidence of unitarity corrections being 
important from the ISR to Fermilab energy regimes. In particular, we sug- 
gest the lack of growth of the single diffraction dissociation cross section [2-5] 
as due to the blackness of central proton-proton and proton-antiproton col- 
lisions. 

In contrast to a very weak growth of the single diffraction dissociation 
cross section in hadronic collisions the energy dependence of virtual pho- 
ton diffraction dissociation appears to be significantly strongc[6-7] than that 
expected from soft physics. In Sec. 3, we argue that this may be due to 
blackness for the soft components of the virtual photon's wavefunction and a 
subsequent dominance of the process by semihard components as suggested 
recently by Gotsman, Levin and Maor[8]. If this is indeed the case it means 
that for the first time one has evidence of unitarity (nonlinear) effects ex- 
tending into the semihard regime of QCD. 

In Sec. 4, wc remark that the new D0data[9,lO] on rapidity gaps between 
jets showing that the gap fraction of events decreases with energy may be 
due to the same physics which slows the growth of the single diffraction cross 
section. If this decrease is indeed due to an energy dependent (decreasing) 
survival probability a similar behavior would be expected for comparable 
photoproduction data involving resolved photons, but such a decrease would 
not be seen in deep inelastic scattering. 

In Sec. 5, we review BFKL searches performed at HI, ZEUS and D0[9,ll- 
13]. Each analysis, using two- jet inclusive measurements at Fermilab and a 
forward jet measurement at HERA, finds some evidence for BFKL behavior 
through an energy dependence which seems stronger than expected from 
leading and next-to-leading order perturbation theory. However, definitive 
results have not yet been achieved. 

In Sec. 6, progress in calculating the next-to-leading corrections in BFKL 
evolution is reviewed[14]. We may be near a rather complete understanding 
of these corrections. A preliminary estimate suggests a substantial reduction 
of the BFKL pomeron intercept. 

Sec. 7 is devoted to a brief discussion of some topics involving nuclear 
reactions. Parametrizations of diffraction dissociation at HERA have been 
successfully used to describe nuclear shadowing in fixed target deep inelastic 
lepton- nucleus scattering[15]. J/ip production in proton- nucleus and nucleus- 
nucleus scattering continues to be an important subject for research. New 
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phcnomenological success in describing all data except for Pb-Pb collisions 
by a simple absorption model[16,17], along with the suggestion that the Pb- 
Pb data may be qualitatively different [18] have made it even more important 
to connect J/%1) production and scattering more firmly with QCD. 



2.1 Total cross sections 

Donnachie and Landshoff[l] have shown that all high energy total cross sec- 
tions for hadron-hadron collisions can be written in the form 



where (Tq depends on the particular hadrons initiating the collision and the 
subleading terms go to zero roughly like (s/so)~^^^- so is an arbitrary scale 
factor while e appears to be universal and of size 



1 -|- e = Qfp is the intercept of the soft pomeron in Regge- language. HERA 
data shows that (2.1) is also true for real photon-proton collisions. Of course 
a growth in energy as fast as that indicated in (2.1) cannot persist at arbi- 
trarily high energies because of limitations required by the Froissart bound 
which does not permit total hadronic cross sections to rise faster than iri^s/so 
at asymptotic energies. The fact that the behavior indicated in (2.1) persists 
up to the highest Fermilab energy region might seem to indicate that uni- 
tarity constraints, which arc responsible for the Froissart bound, are not yet 
effective in the presently available energy region. However, as observed long 
ago[19], this is not the case. If one writes proton-proton total, inelastic and 
elastic cross sections in terms of the S-matrix at a given impact parameter 
of the coUision, S{b), as 



2 Clearly soft 





eRi 0.1. 
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atot = 2 j d%[l - S{h)l (2.5) 

then S{h) depends on b roughly as indicated in Fig.l. (We take S(b) to be 
real for simplicity.) 

For small values of impact parameter S{h) is near zero for proton-proton 
collisions already in the ISR and Fermilab fixed target energy regime. For 
proton-antiproton collisions S{h) is quite small for h < Ifm in the Fermilab 
energy regime. S(b) near zero is a signal that unitarity corrections are large 
though they are not so easy to see in the total cross section because the radius 
of interaction is expanding and the growth of atot is mainly coming from that 
expansion. Monte Carlo simulations [20] of the dipole formulation [2 1,22] of 
the Balitsky, Fadin, Kuraev, Lipatov[23,24] equation for the academic case 
of heavy onium-heavy onium scattering show a similar phenomenon. For 
rapidities less than about 15 the BFKL equation is reliable for the total 
onium-onium cross section, however, for small impact parameter collisions 
important unitarity corrections are visible for rapidities of 6 units. 



2.2 Diffraction dissociation in hadron-hadron scatter- 
ing 

Single diffraction dissociation, illustrated in Fig.2, is given by 

is terms of soft pomeron exchange. Integrating (2.6) over t and over xp < 
0.05, but excluding the proton state = Mp, one gets a single diffractive 
cross section asn- From the Regge formalism one expects asD to grow with 
s as s^*^, but this is not seen in the data as illustrated in Fig. 3, which is 
a simplified version of the more complete plot in Ref.5 where detailed data 
points are shown. A factor of 2 is included in Fig. 3 to account for single 
diffractive dissociation of either of the colliding protons (antiproton). At 
Fermilab collider energies there is a discrepancy of an order of magnitude 
between the Regge fit and the data. 

It seems clear that this discrepancy and the slow growth of ctsd with 
energy signal a breakdown of the Regge analysis when ^/s > 20GeV. I think 
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this breakdown can be expressed in various equivalent ways, (i) A low ener- 
gies inelastic collisions induce, through unitarity, both elastic scattering and 
diffractive dissociation. However, as S{b) goes to zero for small and moderate 
b at high energies, these black regions of impact parameter space only induce 
elastic scattering and not diffraction dissociation. Thus as one increases en- 
ergy the elastic cross section grows rapidly while the diffraction dissociation 
cross section, coming from those regions in impact parameter space where 
S{b) is neither too close to zero or too close to one, grows very slowly, (ii) In 
the Regge language one must include multiple pomeron exchange in addition 
to the single pomeron exchange which is valid at lower energy. This multi- 
ple pomeron exchange gives absorptive (virtual) corrections which slow the 
growth coming from single pomeron exchange, (iii) The "gap survival" prob- 
abihty[25,26] decreases with energy compensating the growth due to single 
pomeron exchange. Although gap survival probability is a concept usually 
used for hard collisions I think the same idea applies to single diffraction 
dissociation, at least in a heuristic way, in hadron-hadron collisions. It is 
likely that (i), (ii) and (iii) are just different ways of saying the same thing. 

3 Deep inelastic lepton-proton scattering 
analogs of the soft physics results 

3.1 An "elastic" scattering amplitude 

Recently, there has been an interesting suggestion as to how to test uni- 
tarity limits in deep inelastic scattering[27]. Of course there is no Froissart 
bound for virtual photon-proton scattering, nevertheless, we have become 
used to viewing the small-x structure function in terms of a high energy 
quark- ant iquark pair (possibly accompanied by gluons) impinging on the 
target proton. Although the quark-antiquark pair is not on-shcll the time 
evolution of the pair as it passes through the nucleon should be constrained 
by unitarity in much the same way that a quark-antiquark pair coming from, 
say, a pion state would be. 

To be more specific, view deep inelastic scattering in the rest system of 
the proton and in the aligned jet (naive parton) model[28,29]. At small x 
the virtual photon, 7*(g), breaks up into a quark and antiquark pair long 
before reaching the proton. The relative transverse momentum of the quark 
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and antiquark is small, of hadronic size jj, ^ 350MeV, while the longitudi- 
nal momenta are Qz and Qz ■ ^ respectively. Because the relative transverse 
momentum is small the transverse coordinate separation of the quark and 
antiquark can be expected to be on the order of a fermi, and the resulting 
cross section with the proton should be of hadronic size. The smallness of 
the overall cross section comes from the small probability, of order fi^/Q'^, to 
find such an aligned jet configuration in the wavefunction of 7*. (More prob- 
able configurations in the 7* have smaller interaction probabilities. While 
the aligned jet model cannot be expected to be a precise model of deeply 
inelastic scattering it should reasonably characterize a significant portion of 
deep inelastic events.) 

The inelastic reaction of this, longitudinally asymmetric, quark-antiquark 
pair with the proton should produce a shadow quark-antiquark pair in the 
final state. If the center of the proton is relatively black to the incoming 
quark-antiquark pair the shadow may be rather strong, as in the hadronic 
case discussed above, and unitarity limits may aleady be reached at present 
energies. The outgoing quark-antiquark pair should show up as a diffractively 
produced state, of mass ~ Q, following the direction of the 7*. Assuming 
that the scattering amplitude of the quark-antiquark pair with the proton is 
imaginary one may reconstruct this amplitude, dropping an i, as 



F(x,b)=/iVey^^ (3.1) 

with b the impact parameter of the collision and p the momentum transfer 
to the recoil proton. ^^§^ is the single diffractive cross section for ~ Q. 

In the present circumstance we do not have good control of the magnitude 
of F near b = 0. However, if the proton is black for central collisions one 
can expect F{x, 0) to show little x-dependence. (Here x plays the role that 
s does for the hadronic collisions discussed above.) The authors of Ref.27 
suggest looking at the b-dependence of 

dinFjx, b) 

- dinl/x ■ ^^-^^ 

Unitarity constraints can be expected to show up as smaller values of Ag// 
near 6 = 0. More quantitatively, unitarity limits at 6 = would give 
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Ae//(6 = 0)<2(«p-l). (3.3) 
This is a clever idea, and it will be interesting to see what the data give. 



3.2 Large mass 7* diffractive dissociation 

n 

The traditional picture of large mass diffraction dissociation at small val- 
ues of X is shown in Fig. 4, where Dokshitzer, Gribov, Lipatov, AltareUi, Parisi 
(DGLAP) [30-32] evolution takes one from the hard scale Q to the soft scale 
IjL where a soft diffractive scattering, represented by soft pomeron exchange, 
occurs. In Fig. 4, one assumes the DGLAP ordering 

/x^ fti « • • • « kl « kl « Q^. (3.4) 

However, this is a subtle process and it is worthwhile looking carefully at 
the argumentation that leads to the size of k^ at the lower end of the DGLAP 
evolution[8, 33-35]. It is convenient to view that evolution proceeding from 
the hard scale Q toward softer scales, a direction opposite to that which is 
usually taken. In Fig. 5, we illustrate the process in two steps: The left-hand 
part of the figure shows the virtual photon wavefunction in terms of its quark 
and gluon components. As in Fig. 4, k is supposed to be the softest gluon and 
Ax± — 2/k± gives the transverse size of the 7* state. The right-hand part of 
the figure gives the diffractive scattering part of the process proceeding by 
gluon exchange from the proton interacting with the octet dipole consisting 
of the gluon k and the remainder of the 7* state. Schematically, one may 
write 



dasD = flux dPr{k±)[l - S{Ax± = 2/k±,h,Y = inl/x)]^(fbdxp (3.5) 
where 

dPr{k±) = ^ (3.6) 

is the probability that the lowest transverse momentum gluon have momen- 
tum k±. Eq.(3.6) shows that gluons with small k± have a small probability in 
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the 7* wavefunction, analogous to what we found earher for low momentum 
quarks in the aligned jet model. 1 — S represents the amplitude for a gluon 
having k±, along with the remainder of the 7* wavefunction, to scatter elas- 
tically on the proton, b is the impact parameter of the overall collision while 
Y = inl/x is the rapidity of the softest gluon with respect to the proton. In 
lowest order, two-gluon exchange, 



xpG{xp, k^) 



^-S^ (3.7) 
^± 

when is large and where an integration has been performed over impact 
parameter, b. Using (3.6) and (3.7) in (3.5) one sees, dimensionally, that /c^. 
cannot be large and this is the logic that has led theorists to take kj_ ^ /i 
and use soft pomeron exchange for the scattering amplitude, 1 — S. 

However, if S is near zero for k± — /i and for b = 0, and this is not 
unreasonable since the S matrix is near zero for small impact parameter 
hadron-hadron collisions, then it is apparent from (3.5) and (3.6) that values 
of k± significantly larger than will be important. Indeed, the values of k± 
that will dominate large mass single diffractive production are those values 
where S is near, but not too close to, one. This is the case since the prob- 
ability in the 7* wavefunction is located in large k±— values. The situation 
here is quite different than for hadron-hadron collisions. In hadron-hadron 
collisions the wavefunction of the incoming hadron is, except for a very small 
part, in the soft physics region. If the (S-matrix is near zero for central colli- 
sions then the inelastic reaction will feed into elastic scattering as a shadow. 
In deep inelastic scattering at small-x when S becomes black there will cer- 
tainly be a similar phenomenon which occurs, and which has been described 
in Sec. 3.1, but, in addition, blackness in the small k± region will allow higher 
values of k± to become effective thus making the process semihard. 

If central impact parameter collisions of 7*-proton collisions are indeed 
black for k± fj, then we would expect the x-dependence of the single diffrac- 
tive cross section, xp-^, to vary more strongly with x than suggested by 
the soft pomeron. If one writes 

xp-^ oc (3.8) 
axp 

then both ZEUS [6] and HI [7] now suggest that n ^ 0.4 rather than the 
n = 2{ap — 1) pa 0.2 predicted by the soft pomeron. If the ZEUS and HI 
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measurements hold up, and n really is near 0.4 in the small (3 region, then 
I think it becomes clear that semihard physics is dominating the physics of 
large mass diffraction dissociation. In that case it is interesting to reexamine 
the "elastic" scattering analyses we described in Sec. 3.1 to see if the proposed 
procedure to measure blackness is also destroyed by the dominance of gluons 
and quarks at higher k^- values. Finally, it should be pointed out that there 
are already rather detailed calculations of the phenomenon, at least for qq 
and qqg components of the 7* wavefunction, which arrived at a value n ~ 0.5, 
not too far from experiment [8] . Have we, for the first time, actually seen the 
long sought after evidence for nonlinearity (unitarity limits) in the semihard 
region of deep inelastic scattering? 

Before leaving this section, it may be useful to again contrast hadron- 
hadron scattering with 7*-proton scattering. In the purely hadronic case 
the energy dependence of the single diffraction cross section is much weaker 
than that predicted by the soft pomeron. We have interpreted this as due to 
blackness in central proton-proton collisions which enhances the elastic cross 
section but suppresses difi^ractive excitation. In 7*-proton scattering, on the 
other hand, the energy dependence (x-dependence) is much stronger than 
that predicted by the soft pomeron. We have interpreted that also as due to 
blackness of the soft components of the 7* now leading to an enhanced role 
for the harder components of the 7*-wavefunction and a resulting stronger 
energy dependence of the cross section. 



4 Rapidity gaps between jets at Fermilab and 



Suppose one measures two jets having comparable but opposite transverse 
momentum along with the requirement that there be a rapidity gap between 
two jets. One might hope that this would be a good process to measure the 
hard (BFKL) pomeron as illustatcd in Fig. 6 [36]. There are, however, at least 
two worries with using this process to measure the hard pomeron. (i) The 
pomeron contribution to the hard quark-antiquark scattering is[37]. 



HERA 



'dt 



4 TT^ exp[2{ap - l)AY] 
4i2 [|aAr,C(3)Ar]3 



(4.1) 
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with Ay the rapidity between the two jets. Here ap — 1 = ^^en2 is the 
BFKL pomeron intercept. The presence of the factor (AY)^ in the denomi- 
nator in (4.1) strongly reduces the effective growth of the cross section with 
Ay making the emergence of the hard pomeron more difficult at moderate 
values of Ay. (ii) Perhaps more serious yet is the fact that the cross section 
for producing two jets with a gap between them depends on the absence of 
a soft interaction between the spectator parts of the proton and antiproton, 
the so-called gap survival probabihty[36]. This lack of factorization makes it 
difficult to make a precise comparison between theory and experiment. 

There is new data from D0[9,1O] and an interesting new analysis com- 
paring the 1800 GeV data with that at 630 GeV. If fgap is the fraction of all 
two-jet events (separated by a given rapidity) with a gap between them then 
D0finds that 

/.ap(630) _2.6±o.6,,„, (4.2) 



/,ap(1800) 

for Ay > 3.8. Thus the gap fraction decreases with increasing energy. 
While this number cannot be directly compared to BFKL dynamics because 
Ay has been taken to be the same at the two energies, while a BFKL test 
should have Ay(1800) - Ay(630) = in^, it does suggest that the survival 
probability has a rather strong energy dependence making BFKL tests more 
difficult in rapidity gap events. It will be interesting to see whether models 
of the gap survival probability can easily accomodate the energy dependence 
in Eq.(4.2)[38]. 

At Fermilab the gap fraction is typically 0.01 while at HERA more like 
0.07. The gap survival probability is much larger at HERA as is natural for 
a point- like 7*. It would be interesting to have a HERA analysis similar to 
that of D0 to see if the energy dependence of the gap fraction is weaker, 
closer to x- independent, than at Fermilab. With respect to the D0data the 
energy dependence of the gap fraction may be reflecting exactly the same 
phenomenon as observed in the energy dependence of the single diffraction 
cross section discussed in Sec. 2. 2. While the inclusive two-jet cross section 
increases at higher energies, because of the growth in the parton densities, 
the energy dependence of the gap cross section is likely to be much weaker 
because of the increasing blackness of central proton- antiproton collisions as 
already seen in the single diffr active cross section. 
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5 BFKL searches 



The hard (BFKL) pomeron or, equivalently, BFKL evolution shows up sim- 
ply only in single transverse momentum hard scale processes. Thus in hadron- 
hadron collisions or in deep inelastic scattering where a soft scale, the size 
of the hadron (proton), is present a special class of events must be taken 
in order to isolate BFKL dynamics. Since this is generally very difficult to 
do experimentally it is perhaps useful to remind the reader why BFKL dy- 
namics is so interesting for QCD and why it is worth the considerable effort 
necessary to uncover it. 

There are at least two important reasons why hard single scale high energy 
scattering is interesting, (i) It is a high energy scattering problem that may be 
soluble, or nearly soluble, (ii) BFKL evolution leads to high parton densities 
and thus into a new domain of nonperturbative, but weak coupling, QCD. 
As parton distributions evolve from a momentum fraction Xi to a smaller 
momentum fraction X2, all at a fixed transverse momentum scale, BFKL 
dynamics gives the rate of increase of those (mainly gluon) densities. This 
evolution is illustrated in Fig. 7. When gluon densities reach a density such 
that on the order of l/o; gluons overlap, perturbation theory breaks down 
and one enters a new regime of strong field, F^,^ ~ 1/g, QCD. While it is 
unlikely that one can reach such densities at Fermilab or HERA at truly hard 
transverse momentum scales one should at least be able to see the approach 
to these high densities through BFKL evolution. 

Inclusive two-jet cross sections at Fermilab and forward single jet inclusive 
cross sections at HERA can be used to measure the BFKL intercept [39- 
43]. These processes are illustrated in Figs. 8 and 9 respectively where ki 
and k2 represent measured jets. In proton- antiproton collisions one chooses 
ki±, kk2± > ^ fixed hard scale, while in deep inelastic scattering ki± is 
chosen to be on the order of Q, the photon virtuahty. For the hadron-hadron 
case 

g(ap-l)Ay 

(^2- jet ^ f{Xl,X2, M^) (5.1) 

while for deep inelastic scattering 
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with Xi and X2 being the longitudinal momentum fractions of the measured 
jets, ap — 1 — ^^^£n2 and the /'s in (5.1) and (5.2) are known in terms 
of the quark and gluon distributions of the proton and antiproton. In (5.1) 
Ay is the rapidity difference between the two measured jets. One can get 
a measurement of ap — 1 in (5.1) by varying AY with Xi,X2 and ]VP fixed, 
and this can be done at Fermilab by comparing the inclusive two-jet cross 
section at different incident energies. In (5.2) one can measure ctp — 1 by 
varying x for fixed Xi and Q^. 

Sometime ago HI [11, 12] presented an analysis showing a jet increasing by 
about a factor of four as x goes from about 3x10^'^ to about 7x10"'^ for 
ki± > 3.5GeV. This is a growth quite a bit faster than given in conventional 
Monte Carlos and much faster than the growth from single gluon exchange 
between the measured jet and the quark- antiquark pair coming from the 
virtual photon. The growth is comparable to that given in (5.2), for ap — 
1 ~ 1/2, however, the comparison is not completely convincing because a 
comparison of partonic energy dependences, from (5.2), with hadron final 
states is not very reliable when ki± is as small as in the HI analyses. 

Recently ZEUS[13] has completed an analysis of this process. Since the 
ARIADNE Monte Carlo gives a good fit to the ZEUS data this Monte Carlo 
is used to unfold the hadronization and thus get a better comparison with 
BFKL evolution. The data agree much better with BFKL evolution than 
with the Born term or with next-to-leading order QCD calculations. A defini- 
tive comparison with BFKL dynamics is hindered by the lack of ability to 
include hadronization corrections along with the BFKL evolution. One can 
hope that the situation will soon improve in this regard. 

A new D0[9] comparing 1800 GeV and 630 GeV data for ki±,k2± > 
20GeV gives ap = 1.35 ± OM{stat) ± 0.22 (syst) when (5.1) is used to fit 
the data. The strength of the D0analysis is that k± > 20GeV which makes 
uncertainties due to jet definition minimal. Weaknesses of the analysis are 
the large systematic error and the smallness of AY, equal to 2, at the lower 
energy. We can hope that the systematic errors will come down in the near 
future. 

Overall, I think the BFKL searches are encouraging but not yet definitive. 
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The fact that all three analyses suggest a strong increase with energy of 
reliable quantities for isolating BFKL effects is certainly positive. An attempt 
will also be made to measure ap — 1 at LEP[44] in the next year by measuring 
the 7* — 7* total cross section. This is a very clean process, although the 
cross section is rather small. 



6 Higher order corrections to BFKL evolu- 
tion 

In general in QCD next-to-leading corrections are very important. It is only 
after next-to-leading corrections have been calculated that scales have a real 
meaning and normalizations can be trusted. In the case of BFKL evolution 
the next-to- leading corrections are also important to show that, in principle, 
corrections to the BFKL answer can be calculated, thus making single scale 
high energy scattering systematically calculable in QCD. 

There has been a long program[14, 45-47], led by the work of V. Fadin and 
L. Lipatov, to calculate the next-to-leading corrections to BFKL evolution 
and it now appears that program may be coming to completion. When the 
work is finished one should get the next correction to ap as well as next-to- 
leading resummations for anomalous dimension and coefficient functions. If 
one writes 



(Xp 



in2a{Q)[l - ca{Q)] 



(6.1) 



then there is the suggestion c may be near 3, a very large correction, although 
there is some work yet to be done before one can accept this number with 
confidence[14]. 

For the anomalous dimension matrix one writes 



7n = 5^7: 



(0) 
ni 



nin — 1) 



+ a^7 



(1) 

ni 



7r(n — 1) 



+ 



(6.2) 



where the first series represents the leading order (BFKL) answer. We should 
soon know the second series, the constants 7^^'* along with similar terms for 
the coefficient functions. When the BFKL corrections are known at next- 
to-leading order we should reap several benefits, (i) A better understanding 
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of the importance of BFKL (resummation) effects in z/W^2 should be possi- 
ble. Recall, that as a two-scale process BFKL dynamics does not directly 
govern the small-x behavior of i'W2. However, BFKL effects are certainly 
present and can be systematically included through resummations such as 
the one indicated in (6.2). When the next-to-leading corrections are known 
we should begin to get a reliable indication of the importance of these resum- 
mation effects in the HERA regime, (ii) The next-to-leading resummations 
should help us to better understand where the operator expansion is valid 
in small-x physics, that is at what x and are coefficient and anoma- 
lous dimension functions sufficiently safe from diffusion effects to be reliably 
calculated perturbatively [47,48] . 

7 Nuclear reactions 

7.1 Nuclear shadowing in deep inelastic scattering 

Nuclear shadowing in deep inelastic scattering is known to be a leading twist 
phenomenon [29] and thus dominated by soft physics, at least at current x- 
values. In the DGLAP formahsm shadowing effects are put into initial parton 
distributions. However, when shadowing is not too strong it can be calculated 
from diffractive deep inelastic scattering using the Gribov-Glauber formal- 
ism as illustrated in Fig. 10. In that figure the left-hand part represents the 
imaginary part of the forward 7*-scattering amplitude which, by the optical 
theorem, is equal to the 7*-A total cross section or, equivalently, h'W2- The 
first term on the right-hand side of Fig. 10 represents the incoherent scat- 
tering off the A nucleons in the nucleus, the nucleons being labeled by Ni. 
The second term on the right-hand side of the figure represents the double 
scattering term which is dominated by diffractive scattering, off nucleon iVj 
in the amplitude and Nj in the complex conjugate amplitude. So long as 
shadowing correction are not too large it should not be necessary to go be- 
yond the double scattering term. Indeed, the double scattering term can be 
obtained from diffractive data at HERA while triple and higher scattering 
terms would involve the scattering of partonic systems in the nucleus, terms 
which cannot be reliably determined. Using a parametrization which fits the 
HERA diffractive data, a pretty good description of fixed target deep inelas- 
tic scattering of nuclei is obtained with the double scattering term giving a 
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shadowing correction of about the right size[15]. 

7.2 J/il) production in proton-nucleus and nucleus-nucleus 
collisions 

Recently there has been much interest, and excitement, about the NA50 
data[18]. In a nutshell one can summarize the experimental situation as 
follows: (i) All P-A and A-A collisions, except for Pb-Pb, look like J/'4> 
production in proton-proton collisions with absorptive final state interactions 
corresponding to a "J/'?/^" cross section with nucleons of 7mb[16,17]. (ii) Pb- 
Pb central collisions have a J/-^ cross section which is significantly suppressed 
with respect to (i)[18] . 

Kharzeev and Satz[49] have suggested a picture that considers the system 
moving though the nuclei, after the hard collision which produces the cc pair, 
to be a {ccg) color singlet system which, if it suffers no reaction with the 
nuclear medium, turns into a, J/ip after the {ccg) system has passed through 
the material. The {ccg) system is supposed to have a size comparable to 
that of the J/-^, but, because of the fact that it looks hke an octet dipole 
formed from (cc)8, and a gluon a cross section of 7mb is natural. This is 
an interesting picture, however, there are a lot of unanswered questions, (i) 
Where does the gluon in the {ccg) system come from? Does it come from 
the hard scattering or is it part of the gluon distribution of the incident 
hadron or nucleus? If it is the latter does this enhance J/tl: production in 
nuclear collisions because there are so many more spectator gluons at the 
impact parameter of the collision, (ii) How does the gluon know what size 
system to form with the cc since the {ccg) does not interact, due to a slowing 
down of the rate of interaction for high velocity states, while traversing the 
material? Why should the relevant size for the ip and ip' be the same? While 
a very interesting, and successful, phenomenology has developed around this 
picture it is important to determine whether the whole picture is reasonable 
from a QCD point of view. 
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